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ABSTRACT 

vj  ^ 

At  the  Fifth  Ionospheric  Effects  Symposium  1 87^ , a paper  was  presented  and  published  that 

described  an  operational  Air  Force  Air  Weather  Service  (AWS),  near  real-time,  ionospheric  specifica- 
tion model  called  Ionospheric  Conductivity  and  Eleotron  Density  (ICED)  . That  version  of  the  model 
specified  primarily  the  northern  hemisphere  midlatitude  and  auroral  zone  ionospheric  regions  and  did 
not  cover  either  the  low  latitude  or  high  latitude/polar  cap  regions.  We  are  making  extensions  and 
modifications  to  ICED  that  provide  a global  specification  — in  effect  a GLOBAL  ICED.  Major  addi- 
tions include  (1)  specification  of  electron  and  ion  density  profiles  (90  to  1000  km)  in  the  high  la- 
titude trough,  auroral  oval,  and  polar  cap  ionosphere  in  both  northern  and  southern  hemispheres,  (2) 
midlatitude  electron  and  ion  density  profiles  that  include  magnetic  declination  effects,  and  <3) 
low  latitude  ionospheric  specification  that  accounts  for  longitude  control  through  differences  in 
vertical  E X B drift  patterns.  Slnoe  this  model  is  driven  by  real-time  data,  we  describe  the  vari- 
ous input  parameters  and  how  they  are  incorporated. 

\ 

INTRODUCTION 

The  Ionospheric  Conductivity  and  Electron  Density  (ICED)  described  by  Tascione  et  al.  [1988] 
provides  near  real-time  specification  of  the  northern  midlatitude  and  auroral  ionosphere  (20  N - 
80°N  corrected  geomagnetic  latitude)  based  on  climatological  models.  The  model  is  driven  by  two 
geophysical  parameters:'  SSNe»»,  an  effective  sunspot  number  (or  ionospheric  index)  determined  from 
the  Air  Weather  Service  real-time  ionosonde  network  observations  and  Qeff,  an  auroral  activity  igdex 
determined  from  satellite  imagery.  ICED  was  modified  to  include  the  low  latitude  ionosphere  (20  S - 
20  N)  by  including  the  Fully  Analytic  Ionospheric  Model  (FAIM)  of  Anderson  et  al.  [1989],  a parame- 
terized version  of  the  Semi-Empirical  Low-latitude  Ionospheric  Model  (SLIM)  of  Anderson  et  al. 
[1987]  combined  with  a modified  version  of  the  Chiu  [1975]  model.  However,  thi3  model  is  only  accu- 
rate for  the  American  longitudes  where  magnetic  declination  is  relatively  small. 

ICED  is  now  being  modified  to  provide  a truly  global  ionospheric  specification  (90°S  - 90°N) 
driven  by  a variety  of  near  real-time  data  including  digital  ionosonde  data,  total  electron  content 
data,  in  situ  plasma  density,  temperature,  and  composition  measurements,  and  satellite  based  ultra- 
violet airglow  and  auroral  emission  measurements.  As  part  of  this  modification  effort,  ICED  is  be- 
ing changed  from  a primarily  climatological  model  to  a more  physically  based  model.-  GLOBAL  ICED  i3 
divided  into  five  regions:;  Northern  High  Latitude,  Northern  Midlatitude,  Low  Latitude,  Southern 
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Midlatitude,  and  Southern  High  Latitude.  The  boundaries  between  these  regions  are  dynamic  in  the 
sense  that  they  are  adjusted  according  to  current  conditions.  For  example,  the  boundary  between  the 
high  latitude  and  mldlatitude  regions  is  the  equatorward  edge  of  the  subauroral  trough,  which  ex- 
pands and  contracts  with  magnetic  activity.  The  northern  high  latitude  region  model  is  in  the  test- 
ing and  validation  phase,  while  the  other  regions  are  still  in  the  development  phase.  The  following 
sections  describe  the  two  high  latitude  models,  the  two  midlatitude  models,  and  the  low  latitude 
model . 

In  all  regions,  we  have  based  the  new  models  on  parameterlzations  of  physically  based  "first 
principles"  ionospheric  models  rather  than  on  statistical  or  climatological  models.  Although  exist- 
ing physical  models  are  still  inadequate  for  aoourately  specifying  the  state  of  the  ionosphere  based 
solely  on  a few  solar  and  geophysical  parameters  and  require  too  many  computer  resources  to  be  used 
in  real  time,  we  believe  that  they  provide  more  realistic  representations  of  actual  ionospheric 
structures  than  oan  be  obtained  from  statistical  models  that,  by  their  very  nature,  tend  to  wash  out 
both  temporal  and  spatial  structure.  In  addition,  the  physioal  models  give  ion  composition  informa- 
tion unavailable  from  models  based  on  ionosonde  measurements  alone.  The  challenge  is  to  find  repre- 
sentations of  these  physical  models  that  can  he  driven  by  real  time  data  using  available  computer 
resources . 


PARAMETERIZING  THE  PHYSICAL  MODELS 

For  this  version  of  GLOBAL  ICED  we  have  chosen  to  parameterize  the  physioal  models  in  terms  of 
geophysical  environmental  parameters  by  producing  databases  of  ion  density  profiles  on  a latitude, 
longitude  (or  local  time),  and  universal  time  grid  for  a range  of  environmental  parameters.  The 
precise  choice  of  environmental  parameters  and  their  values  depends  on  the  ionospheric  region.  Sea- 
sonal variations  are  most  important  for  the  middle  and  high  latitudes.  Solar  activity  is  an  impor- 
tant parameter  at  all  latitudes.  Geomagnetic  activity  is  mainly  important  at  high  latitudes,  while 
-.he  wind-driven  E-  and  F-reglon  dynamos  are  important  for  the  formation  of  the  equatorial  anomaly. 
Since  the  physical  models  Make  intensive  use  of  computer  resources,  we  generally  produce  databases 
tor  high,  moderate,  and  low  values  of  each  parameter. 

Since  eath  database  for  a given  set  of  geophysical  parameters  must  represent  the  ionosphere 
throughout  the  applicable  region  during  a 24  hour  period  in  universal  time,  they  tend  to  be  un- 
wieldy. Therefore,  the  next  step  is  the  production  of  a semi -analytic  representation  based  on  dis- 
crete orthonormal  functions.  (Discrete  orthogonal  functions  have  their  orthogonality  properties  de- 
fi  led  by  sums  on  discrete  grids  rather  than  by  integrals  over  a continuous  interval.  The  functions 
themselves  may  or  may  not  be  defined  between  the  grid  points.) 

The  first  stage  in  the  process  is  the  approximation  of  altitude  profiles  of  ion  density  using 
Empirical  Orthonormal  Functions  (EOF's).  The  generation  and  use  of  EOF's  is  described  in  Secan  and 
Taacicne  (1584)  or  Kutzbach  (1967)  and  references  thertln.  In  short,  the  EOF's  are  the  eigenvectors 
of  the  covariance  matrices  of  the  ion  density  profiles.  If  the  database  consists  of  N altitude  pro- 
files, f(ZjL),  (i»l,2, . . ,M) , the  elements  of  the  covariance  matrix  are 


'ij 


I /n(*i»/n<* ,). 

n=l 


i.J  - 1,2,  ...M 


(1) 


the  high  latitude  databases,  M = 37  and  N = 5760  (20  latitudes,  24  local  times,  and  12  UT's). 

The  eigenvalues  (A  , a=l,2,...M)  of  a covariance  matrix  are  real  and  non-negative,  and  the  ei- 
genvectors, fSB(Zk>.  are  automatically  orthogonal  and  easily  normalized.  (See,  e.g.,  Hildebrand 
(1965).)  He  order  the  eigenvalues  so  that  Am  2 *n  exact  representation  of  the  nth  altitude 

profile,  /n(z).  is  given  by 

^n(zm>  = "‘WV'  m * 1,2, . . ,M,  (2) 

k=l 

where 

H 

“nk  = 1 /n(zm)gk(zm>  (3) 

m=l 


It  cart  be  shown  that  the  RMS  error  due  to  truncating  the  series  (2)  at  K terms  is  minimized 
when  the  are  the  eigenvectors  of  the  covariance  matrix  — rather  than  any  other  set  of  discretely 
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orthogonal  functions  (Secan  and  Tascione  [1984]  and  references  therein). 

For  the  high  latitude  region,  we  have  found  it  best  to  generate  a separate  set  of  EOF's  for 
each  database  rather  than  one  set  for  the  entire  group  of  databases.  We  have  also  found  that  the 
first  six  (out  of  37)  EOF's  produce  good  approximations  and  that  further  improvement  requires  many 
more  EOF's.  Sinoe  the  high  latitude  ionosphere  has  so  much  spatial  and  temporal  structure,  we  ex- 
pect that  the  other  regions  will  require  no  more,  and  quite  probably  fewer,  EOF’s  per  database.  Ex- 
amples of  EOF's  for  0+,  N0+,  and  Oj*  from  one  of  the  high  latitude  databases  are  shown  in  Figure  1. 


Figure  1.  The  first  six  Emplrioal  Orthonoriral 
Functions  for  (a)  0 , (b)  N0+,  and  (c) 
02  derived  from  the  database  for  winter, 
moderate  magnetic  activity,  moderate 
solar  activity,  and  the  BC  convention 
pattern  (By  > 0). 


The  next  stage  involves  the  approximation  of  the  horizontal  variation  of  the  EOF  coeif '.''tents. 
If  we  were  dealing  with  global  databases,  we  would  use  some  form  of  spherical  harmonic  expansion. 
However,  since  we  are  dealing  with  regions  of  limited  latitude  extent,  there  does  not  appear  to  be 
any  special  advantage  to  tils  approaoh,  so  we  have  chosen  to  attack  each  spatial  variable  in  se- 
quence. 

For  each  universal  time  and  magnetic  latitude  the  coefficients,  a^,  of  the  EOF’s  are  repre- 
sented by  Fourier  series  in  magnetic  local  time  (or  longitude).  (Trigonometric  functions  have  the 
useful  property  that  they  are  orthogonal  n uniform  discrete  grids  as  well  as  on  continuous  inter- 
vals.) In  the  high  latitude  region,  we  used  a nine  term  series  (DC  plus  four  cosine  and  four  sine 
terms) . 
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For  each  universal  time  the  Fourier  coefficients  were  represented  by  linear  combinations  of 
discrete  orthonormal  polynomials  in  magnetic  latitude.  These  polynomials  are  generated  for  a given 
latitude  grid  by  the  algorithm  of  Beckman  119731.  We  found  that  a nine  term  series  gave  a good  rep- 
resentation of  the  latitude  variations. 

For  the  high  latitude  model,  the  low  order  polynomial  coefficients  showed  primarily  diurnal 
variations  in  universal  time,  but  the  higher  order  terms  displayed  considerably  more  structure. 
Consequently,  we  chose  to  tabulate  the  coefficients  in  UT  and  use  polynomial  interpolation  rather 
than  to  use  an  orthogonal  function  expansion.  This  may  not  be  necessary  at  middle  and  lower  lati- 
tudes . 


OPERATIONAL  USE  OF  THE  PARAMETERIZED  MODELS 

Once  obtained,  the  semi-analytic  representations  of  the  physical  models  are  U3ed  in  two  differ- 
ent ways.  First,  when  real  time  ionospherio  data  are  unavailable,  they  provide  "best  guesses"  for 
the  current  state  of  the  ionosphere  based  on  estimates  of  the  current  levels  of  solar  and  geomagnet- 
ic activity.  Second,  when  real  time  data  are  available,  they  are  scaled  to  provide  the  best  (least 
squares)  fit  to  that  data,  and  are  used  to  give  the  altitude  profile  shapes.  The  way  in  which  the 
parameterized  model  is  scaled  and  adjusted  by  real  time  data  depends  on  the  region.  One  advantage 
of  working  with  ion  densities  (rather  than  electron  densities)  is  that  the  separation  between  the  E- 
layer  and  F-layer  can  be  made  in  a natural  way:  the  molecular  ions  provide  the  E-layer  while  the 
atomic  ions  provide  the  F-layer. 


The  Northern  and  Southern  High  Latitude  Regions 

The  high  latitude  ionosphere  is  very  dynamic  and  highly  structured  making  near  real-time  speci- 
fication of  electron  density  profiles  extreme  diffloult.  Fortunately,  tor  many  purposes  it  Is  the 
looation  of  the  structures  that  is  of  greatest,  interest  with  electron  densities  of  leaser  signifi- 
cance. We  treat  the  E-  and  F-layers  Independently.  For  the  F-layer,  we  divide  the  high  latitude 
ionosphere  into  three  sub-regions :•  the  Trough,  the  Auroral  Oval,  and  the  Polar  Cap. 

The  dynamic  nature  of  the  high  latitude  ionosphere  is  a direot  result  of  the  dynamic  nature  of 
the  energy  sources  thnt  drive  it:  magnetospherio  oonveotion  and  particle  precipitation.  Further- 
more, the  relatively  long  time  scale  of  the  F-layer  (hours)  implies  that  the  current  state  of  the 
ionosphere  depends  on  recent  history  as  well  as  current  energy  and  momentum  input.  For  this  reason, 
even  a perfect  physioal  model  would  fail  to  accurately  speolfy  the  current  ionosphere  unless  it  had 
accurate  specifications  of  the  reoent  '.'istory  of  the  energy  and  momentum  input.  We  have  adopted  a 
hybrid  approach  combining  both  empirical  and  physioal  models. 

The  physioal  model  we  acopted  is  the  Utah  State  University  high  latitude  ionospheric  model 
(Schunk.  1988],  This  model  calculates  ion  and  electron  densities  and  temperatures  by  following  flux 
tubes  convectlng  through  a mov)  ig  neutral  atmosphere.  An  analytic  version  [Rich  and  Maynard,  1989] 
of  the  convection  model  of  Hepp  er  and  Maynard  {1987]  was  used  to  specify  the  oonveotion  pattern  for 
various  levels  of  magnetic  act)’  Ity.  Ionization  rates  due  to  partiole  precipitation  were  calculated 
using  the  Strickland  et  al.  [19  6]  eleotron  transport  code  and  the  electron  precipitation  model  of 
Hardy  et  al.  [1987].  Photoionic  tion  rates  were  calculated  using  photon  fluxes  from  Hinteregger  et 
al.  [1977]  and  Heroux  and  Hint* rogger  [1978]  (soaled  with  F,0  7)  and  the  weighted  oross  sections  of 
Torr  et  al.  [1979].  * 

The  USU  model  was  used  to  produce  54  databases  containing  N0+,  and  0+  densities  as  func- 
tions of  altitude,  magnetic  latitude,  magnetic  local  time,  and  universal  time  for  a variety  of  envi- 
ronmental conditions.  The  specific  parameters  used  were 

season:.  equinox:  day  82  (March  23) 

June  solstioe:  day  173  (June  22) 

December  solstioe:  day  357  (December  23) 

solar  activity  (?iq  7):  high  (210),  moderate  (130),  and  low  (70) 

magnetic  activity  (Kp) : high  (6.0),  moderate  (3.5),  and  low  (1.0) 

Interplanetary  magnetic  field:  Bz  < 0 (southward) 

B > 0 (BC  convection  pattern) 

By  < 0 (DE  convection  pattern) 

Because  of  the  uncertainty  in  convection  patterns  for  Bz  > 0 (northward),  we  made  no  calculations 
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for  this  case.  Instead,  we  use  the  low  magnetic  activity  cases  to  model  Bz  northward  situations. 

The  first  set  of  54  databases  were  obtained  for  the  northern  hemisphere.  Because  the  relation- 
ship between  the  magnetic  and  geographic  poles  is  different  in  the  two  hemispheres,  a second  set  of 
54  databases  is  being  generated  for  the  southern  hemispheres.  This  set  will  be  treated  in  the  same 
manner  as  the  first  set.  The  northern  hemisphere  databases  have  been  approximated  semi-analytlcally 
using  EOF's,  Fourier  series,  and  discrete  orthogonal  polynomials  as  described  in  the  previous  sec- 
tion. Comparisons  of  the  semi-analytic  representations  with  the  original  databases  are  shown  in 
Figure  2. 
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The  real  time  specifications  of  the  high  latitude  E-  and  F-layers  are  handled  separately.  This 
is  appropriate  because  of  the  disparity  in  time  scales  associated  with  eaoh  altitude  regime.  The  E- 
layer  can  be  handled  in  two  different  ways.  For  post  event  analysis,  when  computer  resources  are 
abundant,  a first  principles  chemical  equilibrium  model  is  used.  Ion  production  rates  can  be  deter- 
mined from  a combination  of  satellite  based  electron  precipitation  data,  satellite  borne  nadir  view- 
ing UV  and  X-ray  intensity  data,  and  statistical  models  of  preoipltation  patterns.  For  real  time 
analysis,  when  oomputer  resouroes  are  soaroe,  the  seml-analytio  representations  of  the  moleoular  ion 
densities  are  scaled  to  provide  a least  squares  fit  to  the  relevant  real  time  data  (mainly  ionosonde 
measurements  of  f0E).  In  addition,  if  information  on  the  looation  of  the  auroral  oval  is  available 
(e.g.,  from  precipitating  particle  measurements  or  from  visible  or  0V  imagery)  the  latitude  scale  is 
adjusted  accordingly. 

Unlike  the  mid-  and  low-latitude  regions,  the  semi-analytic  representation  of  the  high  latitude 
physical  model  is  not  directly  scaled  to  fit  the  real-time  F-layer  data.  Instead,  simple  seml-em- 
plrloal  models  of  the  three  F-layer  subregions  are  used.  These  are  models  of  f-Fj  only,  so  the  pa- 
rameterized physical  model  is  used  to  provide  the  shapes  of  the  0+  altitude  profiles. 

The  F-layer  trough.  Our  semi-empirical  trough  model  is  based  on  the  studies  of  Whalen  [1989], 
There  appear  to  be  three  different  looal  time  features  that  together  make  up  the  complete  trough. 
In  the  afternoon  sector  from  shortly  after  magnetlo  looal  noon  and  fading  out  between  the  terminator 
and  magnetlo  looal  midnight  is  a stable  feature  olosely  related  to  the  ionospheric  convection  pat- 
tern. Beoauae  of  the  long  time  scale  associated  with  F-layer  phenomena,  this  feature  does  not  re- 
spond to  the  instantaneous  convection  pattern,  but  "averages"  over  conditions  prevailing  during  the 
previous  several  hours.  From  the  end  of  the  afternoon  trough  until  sunrise  in  the  F-layer  the  ob- 
served trough  appears  to  be  the  result  of  superimposing  auroral  density  enhancements  on  top  of  the 
naturally  decaying  nighttime  ionosphere.  However,  the  equatorward  edge  of  the  auroral  density  en- 
hancement does  not  necessarily  coincide  with  the  equatorward  edge  of  the  precipitation  because  of 
convection.  Finally,  a morning  trough  sometimes  develops  between  local  dawn  and  local  noon.  The 


analytic  fits  based  on  EOF’s.  £ 

These  profiles  are  from  the 
winter,  high  magnetic  activi- 
ty, high  solar  activity,  By  -q 

positive  database.  The  fit  3 
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precise  conditions  under  which  this  trough  appears  are  not  clear  at  this  tine. 

tn  our  model,  the  equato rward  edge  of  the  trough  is  continuous  in  local  tine  and  is  described 
as  a distorted  circle  with  the  sane  shape  as  the  convection  boundaries  of  Heppner  and  Maynard 
[1987] . The  radius  of  the  circle  is  treated  as  an  adjustable  parameter.  The  width  and  depth  of 
each  trough  feature  (afternoon,  nighttime,  and  morning)  are  independently  adjustable.  In  the  ab- 
sence of  sufficient  data  to  fix  these  parameters,  the  afternoon  trough  is  assigned  a nominal  depth 
of  2.5  MHz  and  a nominal  width  of  5°  latitude,  while  the  nighttime  and  morning  trough  parameters  are 
set  to  reproduce  the  midlatitude  ionosphere.  The  0+  densities  calculated  by  the  midlatitude  models 
provide  "boundary  conditions,"  and  the  depth  parameters  are  calculated  relative  to  those  densities. 

The  auroral  F-layer.  Because  the  auroral  F-layer  is  controlled  by  several  different  and  com- 
peting processes  (e.g.,  precipitation  arid  convection),  and  because  it  is  of  less  interest  than  the 
auroral  E-layer,  we  use  a very  simple  model.  The  equatorward  edge  is  simply  the  poleward  edge  of 
the  trough,  and  the  poleward  edge  is  simply  the  equatorward  edge  of  the  polar  cap.  A single  adjust- 
able parameter  describes  the  peak  fQFj  relative  to  the  f.Fj  value  at  the  pol»r  cap  boundary.  Al- 
though the  potential  to  adjust  the  location  of  the  peak  fQF2  is  retained,  under  normal  circumstances 
it  is  simply  centered  between  the  two  boundaries. 

The  polar  cap  F-layer.  The  polar  cap  is  the  most  difficult  reg^o  to  model,  and  in  the  absence 
of  a dense  network  of  high  time  resolution  ionosondes,  impossible  to  specify  completely.  Under  B 
south  conditions,  convection  brings  high  density  daytime  plasma  into  the  polar  cap.  Because  convec- 
tion is  rarely  a steady  process,  the  usual  result  is  patches  of  varying  size  that  move  through  the 
polar  region.  Under  Bz  north  conditions,  the  polar  cap  ionosphere  is  more  uniform,  but  density  ir- 
regularities at  various  scales  may  still  ooour.  We  attempt  only  to  provide  a crude  model  of  the 
background  ionosphere  and  to  give  statistical  information  on  the  occurrence  of  patches.  If  possi- 
ble, high  time  resolution  digital  lonogram  data  will  be  used  to  estimate  the  interpatch  (background) 
f0Fj  values,  intrapatoh  tJF2  values,  and  the  size  and  spacing  of  the  patches.  The  background  plasma 
density  will  be  obtained  by  scaling  the  CCIR  coefficients  to  digital  ionosonde  observations  in  the 
polar  cap.  In  the  absence  of  such  measurements,  the  current  effective  sunspot  numbe-  will  be  used. 

We  are  currently  in  the  process  of  testing  and  refining  this  part  of  the  model.  Two  prelimi- 
nary testa  of  the  model  using  actual  data  are  shown  in  Figures  3 and  4.  In  both  cases,  only  ground 
based  ionosonde  data  was  used  in  the  parameter  adjustment  process.  Independent  data  was  used  to 
check  the  accuracy  of  the  adjustment  In  reproducing  ionospheric  features. 


Figure  3.  A contour  plot  of  f0F2 
(MHz)  in  geographic  co- 
ordinates based  on  ion- 
osonde data  taken  at 
1200  UT  on  21  December 
1981.  Also  shown  are 
the  locations  of  the 
afternoon  and  morning 
trough  minima  as  deter- 
mined from  DE-2  in  situ 
density  measurements  by 
Santimay  B«3U  (personal 
communication) . There 
were  a total  of  five 
ionosonde  stations  re- 
porting data  in  the 
high  latitude  region. 
The  irodel  required  four 
iterations  to  converge 
to  the  configuration 
shown  here. 
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Figure  4.  A contour  plot  of  f-E 
(MHz)  at  0000  UT  on  21 
December  1981.  Also 
shown  Is  the  equator- 
ward  boundary  of  the 
auroral  oval  (0.3  erg 
cm-2  s"1)  as  determined 
from  NOAA-6  data  by  H. 
W.  Kroehl  (personal 
communication) . There 
were  a total  of  seven 
ionosonde  stations  l e- 
porting  high  latitude 
information,  of  which 
only  one  reported  E- 
layer  data. 


In  both  of  these  cases,  reasonable  agreement  was  obtained  despite  the  limited  data  set  nsec  to 
drive  the  model.  He  anticipate  eve  i better  accuracy  In  normal  operation  with  a variety  of  near  i-oal 
time  data  available.  After  further  testing  and  refinement,  this  model  will  be  validated  against  an 
extensive  set  of  data  intended  to  resemble  as  closely  as  possible  the  expected  complement  of  opera- 
tional data,- 


The  Northern  and  Southern  Midlatitude  Regions 

The  midlatitude  ionosphere  is  less  dynamic  and  less  structured  than  the  high  latitude  iono- 
sphere, Consequently,  it  is  not  necessary  to  resort  to  semi-empirical  subregion  models.  Further- 
more, convection  is  unimportant  ao  the  IMF  orientation  may  be  ignored.  Otherwise,  the  procedure  & u 
very  similar  to  that  used  for  high  latitudes. 

The  midlatitude  physical  model  is  a combination  of  the  photoeleetron  model  of  Jasperse  (198)], 
a photochemical  equilibrium  E-layer  model,  and  a midlatitude  version  of  the  P-layer  (0+)  model  of 
Anderson  (1973], 

The  semi  analytic  repreaentation  of  the  midlatitude  physical  model  is  being  produced  by  the 
same  method  used  for  the  high  latitude  model.  The  real  time  operation  will  bo  somewhat  different, 
however.  A global  ionospheric  index,  roughly  corresponding  to  the  of  the  original  ICED,  will 
be  determined  from  the  complete  set  of  Air  Weather  Service  digital  ionosoncies.  This  index  yili  be 
used  to  scale  the  global  model.  However,  regional  indices  will  also  be  generated  using  ionosonde 
and  other  data  (TEC,  in  situ  plasma  densities  and  temperatures,  UP  imagery,  etc,).  These  indices 
will  allow  the  ion  denaity  scalings  to  vary  with  spatial  location,  resulting  in  a more  accurate  rep- 
resentation of  the  ionosphere  where  data  is  abundant  without  introducing  extraneous  changes  where 
data  is  absent.  As  with  the  high  latitude  data,  the  profile  shapes  come  from  the  theoretical  Bode!,, 
although  the  topside  shapes  will  be  modified  where  appropriate  data  is  available.  In  particular  in 
situ  temperature,  density,  and  composition  measurements  in  the  topside  F-layer  will  be  used  to  ad- 
just the  theoretical  topside  profiles  to  provide  a more  accurate  representation.. 


357 


The  Low  Latitude  Region 

The  low  Latitude  ionosphere  has  many  of  the  characteristics  of  the  high  latitude  ionosphere: 
it  is  dynamic  and  exhibits  spatial  structure.  As  for  the  polar  cap,  we  do  not  attempt  to  apeoify 
the  details  of  the  small  soale  structure  (bubbles,  pluner,  etc.).  The  lnoorpor 'tion  of  FAIM  into 
ICED  provided  limited  low  latitude  capability,  primarily  in  the  region  of  low  magnetic  dedication. 
FAIM  includes  an  analytic  representation  of  the  low  latitude  0+  model  of  Anderson  [1973,  1981]  using 
ion  drift  data  from  Jioamarca  [Woodman,  1 970 j Fejer  et  al.,  1979]  and  assuming  no  magnetic  dedica- 
tion effects.  Ho  additional  real-time  information  (beyond  3SN?ff  and  Qeff)  was  incorporated  in  the 
model.  We  are  using  the  same  model  (combined  with  our  midlatitude  E-layer  model  ) along  with  addi- 
tional data  and  a more  realistio  magnetic  field  as  the  basis  for  the  GLOBAL  ICED  low  latitude  model. 
As  in  the  high  and  middle  latitude  regions,  the  ability  to  adjust  Itself  on  the  basis  of  near  real 
time  data  will  be  an  Integral  part  of  the  model. 

The  semi-analytic  representation  of  the  physical  model  is  obtained  using  the  sume  methods  as 
for  the  high  and  middle  latitude  regions.  A series  of  databases  is  being  produced  for  a range  of 
environmental  conditions.  The  altitude  profiles  will  be  represented  by  a linear  combination  of 
EOF’*,  and  the  coefficients  will  be  represented  by  Fourier  series  in  looal  time  and  discretely  or- 
thogonal polynomials  in  latitude. 

Tn«  real  time  operation  will  make  use  of  in  situ  eleotron  density  measurements  and  UV  airglow 
images  to  determine  the  looation  of  the  equatorial  anomaly  peaks  (or  the  width  of  the  equatorial 
bulge  vh«r  the  anomaly  is  absent)  to  set  the  scale  of  the  electric  field.  TEC  data  and/or  UV  and 
visible  airglow  measurements  will  be  used  to  set  the  soale  of  the  density.  As  in  midlatitudes,  ad- 
justments will  be  made  regionally  aa  well  as  globally  to  produos  the  beet  least  squares  fit  to  the 
ovailafei* 


THE  INTEGRATED  MODEL:  GLOBAL  ICED 

The  five  regional  models  are  to  be  combined  to  produoe  a single  GLOBAL  ICED.  The  midlatitude 
model  will  be  used  to  define  the  values  of  feF,  At  the  equatorward  edge  of  the  trough,  and  these  are 
fixed  at  “boundary  conditions'*  on  the  high  latitude  P-lsyer  model.  Continuity  of  the  0*  profiles  is 
insured  by  taking  a weighted  average  of  the  mid-latitude  and  high-lstitude  profiles  in  the  trough 
wall.  Continuity  of  the  E-layer  (molecular  Ions)  is  Insured  by  doing  a similar  weighted  average 
ovwr  tbs  same  lstit^S*  range. 

The  boundary  between  the  low  and  alddle  latituda  regions  la  lass  sharply  defined.  Because  the 
aidlatitude  and  low  latitude  models  are  basically  tha  same,  tha  principle  difference  being  the  ab- 
sence of  electric  fields  in  tne  midlatltude  model,  the  two  models  must  give  the  same  results  far 
enough  from  the  magnetic  equator.  The  latitude  at  whloh  the  two  models  agree  will  become  the  bound- 
ary between  tha  two  regions  end  will  move  according  to  actual  conditions  as  deteoted  in  the  real 
time  data.  Any  minor  discrepancy  between  the  two  models  will  be  smoothed  out  Oy  using  s weighted 
average  in  the  region  of  overlap. 

Integration  will  proceed  as  soon  as  all  the  regional  models  have  been  tested,  refined,  end  val- 
idated against  actual  data.  Once  Integration  la  oompleta,  GLOBAL  ICED  will  undergo  a final  process 
of  testing,  refinement,  and  validation  before  being  turned  over  to  the  Air  Weather  Service  for  con- 
version to  operation*!  use.  We  expect  GLOBAL  ICED  to  be  oomplete  by  the  end  of  1991, 


FUTURE  DEVELOPMENT 

V"e  see  GLOBAL  ICED  as  the  first  step  in  the  oontinued  evolution  of  ICS).  As  scientific  under- 
standing of  tha  physical  processes  that  control  the  ionosphere  continues  to  improve,  and  so  comput- 
ers continue  to  become  more  powerful,  ve  expect  the  physloal  basis  of  the  modal  to  become  stronger 
and  the  reliance  on  Bsmi-eapirlosl  models  and  semi-analytioal  paraaeteriaations  to  be  reduced.  The 
ultimate  goal  is  a unified  first  principles  model  running  in  real  time  with  periodic  adjustments 
from  re;,)  Uir  wsew-mentg  to  compensate  for  the  uncertainties  inherent  in  any  attempt  to  specify 
the  global  energy  and  momentum  inputs  to  the  ionospherio  system.  We  also  have  begun  the  process  of 
developing  a limited  forecast  model  whose  evolution  will  parallel  that  of  ICED, 
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